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ABSTRACT 

We present simultaneous multi-color optical photometry using ULTRACAM of the transiting exoplanet 
KIC 12557548 b (also known as KIC 1255 b). This reveals, for the first time, the color dependence of the transit 
depth. Our g’ and z’ transits are similar in shape to the average Kepler short-cadence profile, and constitute 
the highest-quality extant coverage of individual transits. Our Night 1 transit depths are 0.85 ±0.04% in z’\ 
1.00±0.03% in g’; and 1.1 ±0.3% in u . We employ a residual-permutation method to assess the impact of 
correlated noise on the depth difference between the z’ and g’ bands and calculate the significance of the color 
dependence at 3.2cr. The Night 1 depths are consistent with dust extinction as observed in the ISM, but require 
grain sizes comparable to the largest found in the ISM: 0.25-1 /.mi. This provides direct evidence in favor of 
this object being a disrupting low-mass rocky planet, feeding a transiting dust cloud. On the remaining four 
nights of observations the object was in a rare shallow-transit phase. If the grain size in the transiting dust 
cloud changes as the transit depth changes, the extinction efficiency is expected to change in a wavelength- 
and composition-dependent way. Observing a change in the wavelength-dependent transit depth would offer 
an unprecedented opportunity to determine the composition of the disintegrating rocky body KIC 12557548 b. 
We detected four out-of-transit u’ band events consistent with stellar flares. 

Subject headings: planets and satellites: detection — planets and satellites: individual (KIC 12557548 b)— 
stars: individual (KIC 12557548) — planetary systems — dust, extinction 


1. INTRODUCTION 


Rappaport et al. (2012 hereafter R12) discovered a 
peculiar, transit-like signal in the Kepler light curve of 
KIC 12557548 therea fter KIC 1255), a K5-K7 dwarf star 
( |Borucki et al.| [2009[ >. The signal is unlike any exoplan¬ 
etary transit curve previously observed, characterised by a 
short pre-ingress brightening phase, a sharp ingress and an ex¬ 
tended egress, with a period of just 15.7 hr. The transit depths 
range from < 0.15% to 1.3%. R12 proposed that the signal is 
caused by a highly irradiated super-Mercury, surrounded by 
an opaque cloud of dust and metal-rich gas constantly replen¬ 
ished by the embedded planet through a thermal Parker-type 
wind feeding a comet-like obscuring tail. In this scenario, the 
underlying planet has an undetectable transit signal, and the 
observed transits are caused by the surrounding cloud of dust 
only. Alternative explanations, such as a dual giant planet sys¬ 
tem or an eclipsing binary orbiting KIC 1255 in a hierarchical 
triple configuration, were also analysed and ruled out. A sim¬ 
ilar signal has since been detected in the light curve of another 
short period system KOI-2700 (Rappaport et al.|2014|. 

A simple cometary dust model for KIC 1255 D was created 
by Brogi et al. (2012 hereafter B12) who estimated an av¬ 
erage particle siz e in the cloud of ~ 0.1 inn. A pseudo 2D 
model by van Werkhoven et al7[( j2013] hereafter W13) was ap¬ 
plied to individual transits in the Kepler short-cadence data, 
suggesting that deeper transits tend to have most absorption 
in the tail, rather than in the disk-shaped, opaque coma sur¬ 
rounding the planet. The pre-transit brightening, attributed 
to forward scattering from the dust cloud, was analysed by 
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[Budaj| ( |2013| hereafter B13), who estimated the particle ra¬ 
dius of the grains in the head of the tail to be a bout 0.1/tm 
to 1 pm, likely decreasing in size along the tail. Kawahara et 
akj (|2013Jl found a correlation between the transit depth and 
tEe activity of the star, pointing to stellar activity being the 
primary energy source driving the evaporation. Simultaneous 
Kepler and /G-band observations of two transits by Croll et 
al. (2014a, hereafter C14) and time-series optical spectropho¬ 
tometry with the 10.4 m GranTeCan (E. Palle 2015, private 
communication) had insufficient precision to detect color de¬ 
pendence in the transit depths at the ~ 10% level. 

In this Letter we discuss high-cadence multi-color observa¬ 
tions of KIC 1255. Our light curves allow us to directly con¬ 
firm the presence of the dust cloud surrounding the putative 
planet and establish that with further similar data it could be 
possible to determine the composition of the scattering mate¬ 
rial. We also briefly discuss flaring events observed in the u’ 
band and compare the morphology of individual KIC 1255 b 
transits with the Kepler observations. 

2. OBSERVATIONS AND DATA REDUCTION 

We observed KIC 1255 on 14, 25, 27, 29 and 31 July 
2013 (civil, beginni ng of the night dates) using ULTRACAM 
( Dhillon et al.|2007) on the 4.2m William Herschel Telescope. 
Dichroic beamsplitters directed incident light into three cam¬ 
eras, each recording a different band. Each light curve was at 
least 4 hr long, covering ~ 1.5 hr of in-transit and > 2.5 hr of 
out-of-transit time, based on the longest transits of R12. We 
used Sloan ( Fukugita et al. 1 1996] ) broadband filters: z\ g’, u 
on the first two nights to observe the widest possible wave¬ 
length range with high signal-to-noise ratio (S/N) to measure 
the color dependence of the extinction and scattering. The 
transit on the second night was very shallow, approaching 
the limits of detectability. Consequently we switched from 
z’ to i’ in the red arm of ULTRACAM for the remaining three 
nights to boost the S/N, while sacrificing some of the wave- 










































2 


Bochinski et al. 


length coverage. We used 6.6 s exposure times throughout 
(co-adding 6 exposures in the u band) to achieve time resolu¬ 
tion up to an order of magnitude higher than the Kepler short- 
cadence data (58.9 sec, W13), achieving S/N > 200 for z’, V 
& g’ and > 50 for u with duty cycle exceeding 99.6% due to 
the use of frame-transfer CCDs in ULTRACAM. All nights 
were photometric, with seeing on average. Guiding er¬ 
rors, due to flexure in the autoguider relative to ULTRACAM, 
caused small offsets in R.A. with ~ 0.3 ,/ (1 pixel) accumulat¬ 
ing per hour, which we manually corrected. The flux showed 
no correlation with the position of the star on the CCDs. 

Frames were calibrated following the standard procedure 
(Howell 2006). with additional fringing corrections applied 
to the first two nights of z’ data. We performed differen¬ 
tial aperture photometry using the ULTRACAM pipelinep 
and the Improved Differential Photometry method (Fernandez 
"" |2012| ), selecting three comparison stars of 
r brighi 


Fernandez et al. 

Hmilar or greater brightness than the target: KIC 12509691, 


KIC 12509684 and KIC 12557520. There were no available 
comparison stars with similar color to KIC 1255, so we per¬ 
formed a subsequent color term correction by fitting a second- 
order polynomial, symmetrical with respect to the meridian 
crossing time, to the out-of-transit section of the light curve. 
3-cr outliers were omitted during determination of the color 
term correction and 8-cr outliers were removed from the final 
results; outliers were fewer than 0.2% of all data. Figure [I] 
shows the resulting light curves for Nights 1 & 2. 

Our light curves show no prominent systematic effects. All 
g’-band light curves for KIC 1255 are shown in Fig. [2] and 
similar apparently systematic-free light curves were obtained 
in other bands and for our comparison stars. We quantified 
the correlated noise in the comparison star light curves and 
in the out-of-transit KIC 1255 data using the time-averaging 
method ( |Pont, Zucker & Queloz|[20061 l. The 10^-5 minute 

range, i.e. the duration of ingress and eg ress is the most 
releva nt timescale for transit photometry ( |Carter & Winn| 
2009 ); the red noise on this timescale in our KIC 1255 Night 

1 light curve is on the order of 6%, 9%, and 30% of the white 
noise in z’,g’, and u respectively. The noise levels on other 
nights were comparable, and in the comparison stars the red 
noise levels were slightly lower, especially in u . 


3. RESULTS AND ANALYSIS 
3.1. Transit Characterisation 


Following Haswell et al. (2012|) we used a simple Man- 


g Haswel 

del & Agol (2002) planet transit model to determine R p /R 
and the mid-transit time of KIC 1255 b in each band for each 
night. We set the remaining planetary parameters according to 
W13 and used stellar parameters from the Kepler Input Cata¬ 
logue (Bro wn et al.|20lf) . Quadratic limb-darkening coeffi- 
cients were interpolated from Claret & Bloe men|( 2011) using 
the EXOFAST routine by |Eastman et al.|(|2013). We know 
the transits of KIC 1255 b are not attributable to the unchang¬ 
ing circular cross-section of an opaque planet, but this simple 
modelgives quantitative depth and timing information (see 
Tablem), which can be easily reproduced. 

The Mandel & Agol fits to the Night 1 z’ and g’ give resid¬ 
uals which show clear systematic deviations, indicating as 
expected that the dust is distributed differently to a circular 
cross-section. We also fitted scaled versions of the average 

4 http://deneb.astro.warwick.ac.uk/phsaap/software/ultracam/html/ 
index.html 


Kepler short-cadence profile (C14). In this fitting we allowed 
the time/phase offset to vary freely and applied a scale fac¬ 
tor to the flux, so the fractional depth, d, at the deepest part 
of the profile was free to vary, but the duration of the transit 
profile was fixed. Like our g’ and z’ data, the Kepler short- 
cadence profile has some discontinuities in the slope of the 
egress, and provides fits which do not exhibit systematic pat¬ 
terns in the residuals. It is possible that the very smooth egress 
in the phase-folded Kepler long-cadence curve (R12, com¬ 
pared to our data and to the Kepler short-cadence curve) can 
be partially attributed to the smearing effects of the long Ke¬ 
pler exposures. Our light curves, particularly the g’ band, 
which is similar in wavelength coverage to the Kepler band, 
are very similar in shape to the scaled Kepler short-cadence 
profile. Figure HI shows the best fitting scaled Kepler short- 
cadence profilesfor Nights 1 & 2, while Figure [3] shows the 
corresponding confidence intervals for Night 1. Following the 
rationale in Carter & Winn|( j2009| ), the confidence intervals in 
Figure [3] a) and Table [T] include the contribution from the red 
noise. 

We used the (only published) fiducial phase from R12 and 
the (most precise published) period from W13 to calculate 
phase 0.0 independent of our data. Figure HI shows time rela¬ 
tive to this. Where our unbiased search of depth-timing space 
identifies a transit at the appropriate phase, this is sugges¬ 
tive of a genuine detection. We detected appropriately phased 
transits on Nights 1 and 2; the best fit parameters are reported 
in Table [T] The Night 2 g' light curve (Figures []]& [2} clearly 
has a noise level similar to tbe reported transit depth, and the 
correlated noise can conspire to move the out-of-transit level 
and the in-transit level in opposite directions, so while the tim¬ 
ings suggest this is a genuine detection we caution against 
placing too much faith in the reported depths. 

For Night 1 the reduced ~/ * 1 2 * 4 values for our Mandel and Agol 
model fits were 0.872, 0.853, and 0.800 with 4036, 4036 and 
671 degrees of freedom in z’, g’ and li respectively; the cor¬ 
responding values for the scaled Kepler short-cadence profile 
fits were 0.785, 0.695, and 0.785. The scaled Kepler profile 
obviously provides a superior fit to the z’ and g' data. The u 
light curve (Figure [TJ clearly does not strongly constrain the 
shape of the transit. 

We report the mid-transit time for the Mandel and Agol 
fits, while for the Kepler profile fits we follow C14 and re¬ 
port /transits the time of flux minimum. Because the profile is 
skewed, /transit occurs after the mid-transit time; the Night 1 z’ 
and g’ values in Table[I]show offsets of ss 0.003 days, which 
is within the uncertainties in the u’ values. The depth param¬ 
eter, d , is the fractional depth at the deepest part of the fitted 
Kepler profile. For ease of comparison with the Rp/R* param¬ 
eter, we tabulate the values of d^. Again, there is a systematic 
offset: the broader, flatter floor in the Mandel and Agol transit 
profile sits above the minimum of the observed profile, while 
the Kepler profile matches it well. This causes the system¬ 
atic offset between the tabulated Night 1 depth parameters of 
« 0.01 in all three bands. 

On Nights 3, 4, & 5, whose g’ band light curves are shown 
in Fiure. [2] we did not recover appropriately phased tran¬ 
sits. We estimated upper limits on R p /R* on these nights 
by attempting to recover transits using the method described 
above, and quoting the largest value of R p /R* identified: 
0.031 in i\ 0.027 in g’ and 0.087 in u’ (averaged over three 
nights). These ‘transits’ were at the wrong phases, indicating 
the algorithm chose noise features. Any actual transit deeper 
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KIC1255 - Night 1 


KIC1255 - Night 2 



-0.10-0.05 0.00 0.05 0.10 0.15 



-0.10-0.05 0.00 0.05 0.10 0.15 
time (MJD-56487.986) 



-0.10-0.05 0.00 0.05 0.10 0.15 



-0.10-0.05 0.00 0.05 0.10 0.15 
time (MJD-56499.101) 


Figure 1 . Upper panels: z’ (top), g ’ (centre) and u ’ (bottom) light curves of KIC 12557548 from 14 & 25 July 2013; z* and w’ are offset vertically by +0.03 and 
-0.03 respectively. The scaled Kepler short-cadence profile fits are shown in red. four-minute bins (z’ and g ’) and 24 minute bins (m’) are shown in black. Grey 
points ( u ’) show an alternative binning of the same data. Raw data points, shown in orange for z’ and green for g' are omitted for Lower panels: g’/z’ curves. 
Times are relative to a zero-point calculated from the fiducial phase of R12 using the period of W13. 


than these values would have been found if present. 

If the obscuring dust is clumpy, we might expect an in¬ 
creased scatter in flux values during ingress and egress as in¬ 
dividual clumps cross the stellar limb. We searched for this 
effect but found no significant evidence for it. It may, how¬ 
ever, emerge if a larger number of light curves of the quality 
of our Night 1 data are analysed. 

Our primary result (Table ITJ is the color dependence of 
the transit depth: the ratio ofthe square roots of the transit 
depths in g’ and z’ (i.e. the parameter analogous to R p /R *) is 
1.084 (+0.024, —0.029), where the uncertainties were propa¬ 
gated from the best fit model depths. To assess the robustness 
of this, we employed residual permutation statistics (Winn et 
al. 2009) with two-sided 68.27% confidence intervals based 
on the studentised pivotal method ((Carpenter & Bithell|2000|, 
concluding that the ratio of the square roots of transit depths in 
g’ and z’ is 1.084 (+0.039, —0.033). As none of the 807 real¬ 
isations of the residual permutation method produced a depth 
ratio of 1 (Figure [3j b)), we reject the zero color difference 


hypothesis with probability 99.88% or 3.2a. 

To assess whether the transit depth, which can be influenced 
by limb darkening, is a good proxy for light loss due to ex¬ 
tinction in a given band, we compared the g’/z’ depth ratio 
described above to the ratio of the equivalent occulting areas 
A„i/A z i in these two bands. To find A„i & A z i, we modelled 
the occulting dust tail as a cloud with variable transparency 
extending 0.075/?* above and below the orbital plane and 60° 
following the orbit of the planet (B13). The dust density pro¬ 
file of the tail, inversely proportional to its transparency, is: 

P(0)=P( 0)e Ce/ * (1) 

i.e. Equation 4 from B13; where 9 is the angle from the be¬ 
ginning of the arc, p(9) is the density along the arc and C 
is a constant. C = —25 produces a good match to the tran¬ 
sit light curve in all bands, and is consistent with B13. The 
transparency along the tail was adjusted by varying p( 0) in 
each band separately, until the light curves matched the ob¬ 
servations. The square root of the resulting A g i /Aj ratio is 
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Figure 2. Our g’-band light curves. Vertical offsets of (n — 1) X 0.035, 
where n is Night number, have been applied. Dashed line: transit time from 
R12 and W13 (see text). Dotted lines: transit duration (R12). 



1.077 ±0.001 (the uncertainty arises from a reasonable range 
of limb-darkening coefficients), well within the confidence in¬ 
tervals for the square root depth ratio. This establishes that the 
impact of limb darkening on our results is minor in compari¬ 
son to other sources of error. 

The presence of occulted and unocculted starspots can 
modify the transit depth and possibly influence the derived 
dependence of the planetary radius on the wavelength (Pont 


et al. 


2008). To assess this, we looked at analogous multi¬ 
wavelength observations of a well studied mid-K star, HAT- 
P-20 (IBakos et al.|201 1 j ), and used Figure 4 in Ballerini et al. 


(2012) to calculate the effect of starspots on the theoretical 
-20 g’/z’ transit depth ratio. The maximum fractional 
light loss_due to starspots was set at 4% as predicted by |Croll| 
et al. (|2014b]> for KIC 1255. We found that starspots can ac¬ 


count for only up to 0.46% of our observed g’/z’ square root 
depth ratio: insignificant in comparison to other sources of 


3.2. Color Dependence and Dust Modelling 

We observed differences in transit depths between the u’ , 
g’ and z’ filters on Night 1, with an increase in depth toward 
shorter wavelengths (Figures |T] & [3] Table [lj. This effect is 
consistent with extinction from the putative dust cloud sur¬ 
rounding the planet (R12), as we will show quantitatively. De¬ 
tections on Night 2 (Figure [I] Table [T} also show an increase 
in depth towards shorter wavelengths. 

We modelled our Night 1 scaled Kepler depths with the 
|Cardelli, Clayton & Mathisj ([1989] hereafter CCM89) inter¬ 
stellar extinction law to find that the extinction depths are 
broadly consistent with the wavelength-dependence of the 
reddening seen in the ISM (Fig. HV This is strong indepen¬ 
dent and direct evidence in favor or the dust-cloud model for 


the transits. We constrained Ry to the highest values found in 
the ISM: R y sa 5.3. As Ry increases with the mean dust grain 
size ( Mathis & Wallenhorstl|l98T| l, the grain sizes we found 
in the KlC 1255 dust cloud on Night 1 are comparable to the 
largest found in the ISM (0.25-1 /Tin depending on the grain 
composition, Mathis et al.|1977[ ). We normalised and plotted 
the two transits from C14 in the inset to Figure |4j their color 
dependence is consistent with our fit, as well as with a con¬ 
stant transit depth. Note that the C14 transits were observed 
on a different night to our observations and hence a different 
grain size distribution might apply. The transits vary because 
the dust cloud is varying, and it is difficult to imagine a sce¬ 
nario where the irradiated dust cloud changes on planet-sized 
scales while the dust grains comprising it retain an identical 
size distribution. 


3.3. Flares 

We observed two potential flaring events in the u’ band 
on Night 1 (Figure |T|), characterised by a sudden increase 
in brightness by approximately 1.0% and 1.5% with a sim¬ 
ilarly abrupt drop back to base level. On each of Nights 2 
and 4 there were single similar events with amplitudes 1.0% 
and 1.5% respectively. None of our comparison stars showed 
anything similar, suggesting these are genuine astrophysical 
signals originating from KIC 1255. The atypical flare shape 
can be attributed to the faintness of the target coupled with 
the low effective cadence of our u’ band data hiding the expo¬ 
nential decay. On Night 1, the second flare overlaps in phase 
with a slight brightening in z’ band. Similarly on Night 2, a 
bright flare in the u’ band overlaps in phase with a brighten¬ 
ing in both the z’ and g’ bands. To quantify the amount of 
extra emission expected in the u’, g’ and z’ bands, we mod¬ 
elled continuum or ‘white light’ flare emission as a blackbody 
curve with a temperature of 10,000 K and a variable fill fac¬ 
tor, defined as the effective surface area of the star covered 
by the active region ( [Kowalski et ak|2010 i. We added this to 
the model spectral energy distribution for KIC 1255 (|Howarthi] 


Table 1 

Parameters derived from fitting the Mandel & Agol model and the averaged 
Kepler profile. 


Filters 

Night 1 

Night 2 


R P /R, 

d 2 t 

R P /R* 

d \ t 

z’ 

o.o826 ±8:8811 

0 0999 +0 0020 
u.uvzz _ 0- 0016 

0 023 + 0 006 
V.uzj) -0.006 

0 029 + 0 007 
u.uz? _ 0 007 

g’ 

0.0B91 

0 0999 + 00013 
v vyyy - 0.0015 

0 035 + 0 004 
-0.004 

0 037 + 0 006 
u.u^/ _o.o06 

u ’ 

0 090 + 0 016 
u.uvu _ 0 .020 

0-104 Iflojo 

0 097 + 0 031 
u - uy/ -0.043 

0 137 + 0033 
UA:>/ -0.031 

g’/z’ 

1 079 + 0027 
lu/y -0.021 

1-084 IK* 

1 52 +0 - 43 
l - DZ - - 0.43 

1 98 +0-37 

i-ZQ _o.37 


mid-transit 

^transit 

mid-transit 

^transit ^ 

z’ 

56487.9858 

+0.0004 

-0.0005 

56487.9886 

+ 0.0010 

-0.0007 

56499.103 

+ 0.010 

- 0.001 

56499.110 

+0.006 

-0.005 

g’ 

56487.9866 

+0.0004 

-0.0004 

56487.9901 

+0.0004 

-0.0004 

56499.092 

+0.006 

-0.001 

56499.101 

+0.005 

-0.002 

m’ 

56487.986 

+0.007 

-0.004 

56487.992 

+0.004 

-0.007 

56499.106 

+0.007 

-0.004 

56499.113 

+0.005 

-0.005 


^ Square root of the transit depth from Kepler averaged profile fitting. Com¬ 
parable to the Rp/R* value. 

^ Transit time as defined in ( 


Croll et al. 


2014a [I 


Alternative result derived from the residual-permutation method is 


1.084 


+0.039 

-0.033 
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time (MJD - 56487.986) 



g'/z' depth 172 ratio 


Figure 3. (a) Contour maps from our fits of the scaled Kepler short-cadence profile to the Night 1 z* (left), g’ (central) and u ’ (right) bands. Black lines give 

confidence intervals and red dashed lines show la uncertainties, (b) Histogram of 807 residual permutati ons for the g’/z’ d ep ths rat io. The dashed line shows 
the best fit and the dotted lines show 68.27% confidence intervals based on the studentised pivotal method {Carpenter & Bithell|2000) . 


CCM extinction curve 



Figure 4. The CCM89 ISM extinction law fitted to our Night 1 transit 
depths. The fit included two free parameters: E(B-V) and Ry. The best fit 
values were E(B-V)=0.056 (used throughout) and Ry = 5.3 (solid line). Ex¬ 
tinction laws for values of Ry ranging between 2.50 and 9.00 are also shown. 
The inset shows our best fit CCM law and C14’s photometry normalized to 
the same transit depth; it is consistent with our fit. C14’s Night 2 vertical 
error bars (dashed lines) have been offset slightly in wavelength for visibility. 

2011J) and convolved the resulting artificial spectrum with fil¬ 
ter responses for the SDSS and Kepler bands. The best fit to 
the observed ~ 1% n’ band flare signals was achieved for a 
fill factor of 0.013%, with a surplus flare to stellar flux ratio 
in the u\ g\ z' and Kepler bands equal to 1.0%, 0.3%, 0.1% 
and 0.5% respectively. This shows that flares originating on 
KIC 1255 are expected to be brightest in the u band, with 
signal in the remaining SDSS bands comparable to our noise 
level. 

The calculated surplus emission in the Kepler band is con¬ 
sistent with ‘white light’ flares observed on other K dwarfs 
(between 0.1% and 1,0%; |Walkowicz et al.|201 1) . If flares on 
KIC 1255 are similarly ‘white’ in nature, they should be bor¬ 
derline visible in the Kepler short-cadence photometry. We 
inspected all currently accessible short-cadence Kepler pho¬ 


tometry (Q13-Q17) and found several similar brightening ef¬ 
fects with an amplitude of up to 0.5% and duration between 
30 and 90 minutes. These are consistent with our predictions, 
but also comparable to the Kepler photometric precision at 
~ 0.4%. Additionally, a 1200 days long cadence Kepler light 
curve of KIC 1255 b shows at most three significant outliers 
that could be consistent with flaring events or simply outliers, 
the strongest having an amplitude of 3. 5% (S. R appapo rt & B. 
Croll 2015, private communication). Kawahara et al.]( |2013| 
found that KIC 1255 is capable of producing even brighter 
flares, with a brightness increase of up to ~ 6% on timescales 
of 30 minutes. 

4. DISCUSSION 

We obtained multi-color optical light curves of the transits 
of KIC 1255 b. For the first time these transits are measured 
with sufficient precision to discern features of amplitude 0.2% 
and duration ~ 4 minutes in an individual transit light curve. 
This permitted robust quantitative characterisation of the tran¬ 
sit shape and depth on Night 1. We clearly detect a difference 
in depth between the z’ and g' bands: visible in the lower 
panel of Figure [I] and quantified in Table 1 and Figure [3] We 
rule out the effects of limb darkening and starspots as signif¬ 
icant contributions to this depth difference. The color depen¬ 
dence of the transit depth is consistent with extinction due to 
dust with grain sizes 0.25/im < a < 1 ,um, where a is the grain 
radius. 

Eleven nights later, on Night 2 we detected much shallower 
z’ and g' transits, with no detection of color dependence in 
the g' and z’ depths (Table [TV The changing transit depth 
is probably primarily causedby changes in the quantity and 
spatial distribution of dust due to a limit-cycle behavior of 
the dust injection mechanism (R12). Clearly the dust cloud 
produces less extinction during the shallow and undetected 
transits than it did on Night 1. It would be very interesting 
to examine the color dependence of the transits at a variety 
of depths: this will reveal whether the grain size distribution 
changes with transit depth. 

Figure 13 of C14 shows dust grain extinction efficiency as a 
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function of X = 2na/\ for grain size a and wavelength A. For 
X >~ 2 the efficiency plateaus at its peak value irrespective of 
composition. For smaller values of X the extinction efficiency 
drops off in a composition-dependent way. As X <x A -1 , this 
lessening of extinction efficiency with decreasing grain size 
would be most pronounced at long wavelengths. With further 
observations of quality comparable to Night 1, with transits 
of varying depths but deep enough for secure detections, we 
have the prospect of identifying the composition of the dust 
grains by matching the observed wavelength-dependent depth 
changes to one of the curves in Figure 13 of C14. In our 2013 
campaign, we were extremely unlucky to observe rare shallow 
transits with depth < 0.07% (in g’) on four of the five nights 
we observed. Only two such quiescent periods were observed 
by W13 in years 2009 - 2013, amounting to ~ 3% of total 
Kepler observing time. 

We have presented direct evidence that the variable transits 
in KIC 1255 b are due to an evolving cloud of dust attributable 
to mass loss from a disintegrating low-mass rocky planet. 

Analysis of a larger number of transit light curves poten¬ 
tially offers the opportunity to measure the composition of 
the grains. This is an unprecedented opportunity to learn the 
make-up of a rocky extrasolar planet: further observations are 
urgently required. 
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